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3.1 Introduction 

Recent decades have witnessed numerous experimental procedures for creating a very sharp spatial 
temperature inhomogeneity in various materials. Thus, temperature gradients of up to 10 K/cm can be 
established by laser radiation on semiconductors, electron-beam material processing, welding, etc. Large 
electron temperature gradients can also be achieved by heating up semiconductor charge carriers with 
microwave field ' ; hence the relevance of investigation into semiconductor properties under large 
temperature gradients. 

The tendency toward miniaturization of thermoelectric power converters also leads to the creation of 
large temperature gradients in semiconductor materials. Consequently, the task arises to study 
transport phenomena in semiconductors under extreme temperature conditions. This situation will 
arise in the use of submicron and nanosize structures as thermoelectrics. Thus, in a thermoelectric 
generator based on a quantum super lattice, for example, temperature gradients up to 10 5 K/cm are 
established. 3 Under these conditions, an extremely high thermoelectric figure-of-merit was obtained: 
ZT - 2.4 at T = 300 K. 3 

In addition, new effects can be expected under large temperature gradients. One is the Benedicks 
effect, which has been described in a number of experimental and theoretical papers. The origin of 
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FIGURE 3.1 Conditions of the Benedicks effect observation: (a) sample shape, (b) temperature distribution T(x) 
along the sample. 

this effect lies in the thermo-electromotive force (EMF) in a homogeneous isotropic sample where two 
opposite temperature gradients of different values are created — the ends of the sample being maintained 
at equal temperature (Figure 3.1). 

In this chapter the mechanisms of the effect of a temperature gradient on electric current transport 
phenomena in semiconductors are discussed and a novel mechanism for this influence is proposed, 
namely, the heating up of charge carriers under the influence of a large temperature gradient. 

3.2 Influence of a Temperature Gradient on Transport 
Phenomena in Semiconductors 

The influence of temperature gradient on charge carriers transport starts to be manifest when 
characteristic temperature changes become comparable to the mean free path of electrons or to any other 
characteristic lengths. 

The characteristic length Lj of temperature change is equal to 1 ' 

T 

(3.1) 



Ivrl 



where VT is the temperature gradient. 

The following inequalities hold true for the principal microscopic lengths in elastic or quasielastic 
mechanisms of charge carrier scattering 1 ' : 



L r » L e » L p » L ph 



(3.2) 



where L r is the recombination length, L e is the current carrier cooling length, L„ is the momentum mean 
free path of carriers, and L„ h is the thermal phonons mean free path. 

With increasing temperature gradient, the length L T can first of all be comparable to the largest of the 
lengths in inequalities (Equation 3.2), the recombination length 

L r > L T (3.3) 

Tauc 7 theoretically investigated the case of Equation 3.3 when studying the Benedicks effect. According 
to the Tauc theory, under the condition of Equation 3.3, the minor current carriers diffusing into the low- 
temperature region of the sample at a sudden temperature drop are not in equilibrium. This causes a 
change in the semiconductor properties in the region of large temperature gradients, for example, a change 
in the thermo-EMF. However, a change in thermo-EMF due to large temperature gradients in 
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this case proved to be very small and can be completely neglected in materials with a large current carrier 
concentration. 

With further increasing the temperature gradient, the characteristic length L T can become comparable 
to the cooling length L E : 

It < K (3-4) 

In this case, it is natural to expect the appearance of new electron transport phenomena caused by the 
analogy between strong electric fields and large temperature gradients. It is appropriate to introduce the 
concept of specific heating of the electron gas in semiconductors under large temperature gradients. 
Thus, for a homogeneous isotropic semiconductor with a standard band structure the expression for 
current density in the presence of electric field and temperature gradient is given by the kinetic integrals: 



■in 



3jirmT 5/2 



\e^ s 3 \e° ,T de " Vr[ | JJ e 3 V E/T ds - J \] s 5/ V /r d £ ] j (3-5) 



where n, m, e, and r p are the concentration, the effective mass, the energy, and the electron momentum 
relaxation time, respectively, and £is the chemical potential (for simplification offormulas in Equation 3.5 
and in all subsequent equations of Chapter 3, we accept the system of units in which the Boltzmann 
constant K B — 1 and the electron charge e — 1). 

Provided eddy currents are not present in the case under consideration, in an open circuit an electric 
field is created: 



(l"" 1 "^ 



VT (3.6) 



where t„ ~ e" 1 . Using Equations 3.5 and 3.6 the specific density of partial current (electron current with 
the energy from e to [e + de]) can be written: 

1T = 3^1*° lr""'"2J Vr (3 - 7) 

It follows from Equation 3.7 that even at zero macroscopic current the partial currents are nonzero, 
because the difference in the parentheses goes to zero with the only fixed value of e. As long as the partial 
currents are temperature dependent, and the temperature depends on coordinate x, the number of 
e-electrons entering the element of the sample [x, x + cbc] is not equal to the number of e-electrons 
leaving the element. Thus, the number of e-electrons in this element of the sample will change. 
Consequently, the symmetric part of the distribution function (concerning to the quasimomentum) will 
also change. The change in the distribution function can be interpreted as specific electrons heat up 
under large temperature gradients. Naturally, the heating up will be significant if the temperature 
gradient is large enough. 

If the conditions for establishing the electron and phonon temperatures are provided, the nature of the 
electrons' heat up under large temperature gradients can be also understood from macroscopic 
considerations. ' 

As long as L E continues to be the second largest microscopic length in semiconductors next to the 
recombination length L r , the investigation of phenomena caused by the heat up of electron gas under 
large temperature gradients will be of primary practical interest. 

The cooling length is also a very important parameter for studying the heat up of electrons; cooling 
lengths at room temperature for a number of materials are listed in Table 3.1. ' Cooling lengths, as a rule, 
increase with a decrease in temperature. 
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TABLE 3.1 Cooling Lengths for Different 
Materials 



Material L e (cm) 

Ge 1.7 X 1(T 4 

Si 5.9 X 1(T 5 

GaAs 4.0 X 1(T 4 

InSb 1.0 X 1(T 4 

BiTe 5.5 X 10~ 5 

CdS 5.6 X 10~ 6 

CdSe 2.7X10~ 5 

HgSe 9.1 X 10~ 4 

GaP 7.1 X 10~ 6 

InP 2.3 X 10~ 4 

PbTe 7.7 X 10~ 5 

HgTe 1.0 XHT 3 

Sn 4.0X10~ 6 

Pb 2.6 X 10~ 6 

In 5.6X10~ 6 

Cu 2.0 X 10~ 6 

Ag 2.1 X 10~ 5 

Au 1.2X10~ 5 



3.3 Calculation of the Influence of a Large Temperature 
Gradient on the Properties of Semiconductor Materials 

Consider a homogeneous, isotropic nondegenerate semiconductor with a standard band structure. To 
describe the influence of a large temperature gradient, it is necessary to take into account the self- 
consistent behavior of electron and phonon subsystems and the electric field in the sample. Consider two 
variants of the problem formulation. 

3.3.1 High Current Carrier Concentrations 

It is assumed that the current carrier concentration is sufficiently high; i.e., for the current carrier 
subsystem the electron temperature approximation can be used. The respective inequalities are written in 
Ref. [8]. Under these assumptions the equations of electron and phonon energy balance can be derived 
from the original kinetic equations by respective integration ' : 

div q e =]E* -P(T e -T p ) (3.8) 

div q p = P(T e - T p ) (3.9) 

where E" = — Vip, cp is the electrochemical potential, T e and T p are the electron and phonon temperatures, 
respectively, and the function P — P(T e , T p ) is responsible for electron-phonon energy interaction. The 
obvious form of P function for various scattering mechanisms can be found in Ref. [8]. For electron and 
phonon densities of heat flow q ep and for current density j the phenomenological relations can be used: 

i=o<r e )£*-a<r e ) a (r e yvT e (3.10) 

q e = -K e (T e )VT e + n(T e j] (3.11) 

q e =-K p (r p )vr p (3.12) 
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where cr(T e ), a(T e ), K e (X e ), 77(X e ) = T e a(T e ) are the electric conductivity, the thermoelectric coefficient, 
the electron thermal conductivity, respectively, and the Peltier coefficient, kJT„), is the phonon thermal 
conductivity. Equations 3.10 to 3.12 should be supplemented with the current continuity condition in 
Equations 3.8 and 3.9 

divj = (3.13) 

and the system of Equations 3.8 to 3.13 with boundary conditions. 

Experimental results ' prove convincingly that the thermal boundary conditions for electrons and 
phonons can be assigned separately. For this reason sufficiently general thermal boundary conditions 
have been assumed (in a one-dimensional case): 

&,pl*=-a ~ -T?eTp( T e,p ~ ^l)!^-,,, q e J x=a = 17e,p(X eiP - T 2 )\ x=a (3.14) 

where 17^ are the phenomenological coefficients of heat exchange depending on conditions in 
contact areas. 

This written system of nonlinear differential equations shown above in the general case cannot be 
solved analytically. Its analytical solution cannot be obtained even in a single-temperature approximation 
(T e = T p ) for the simplest boundary conditions. 

Assume that the relative temperature drop at the ends of the sample is small compared to the average 
temperature: 

T 2 ~ T, T 2 + T, 

y= r « 1. T* = -^-^ (3.15) 

Here, T 12 are determined by the respective boundary conditions (3.14). Under such assumptions the 
problem can be divided into two parts: 

(A) Determination of the electron and phonon temperatures and electric fields in the sample volume. 
For this purpose the following solution methods were developed and used: the modified method of 
perturbation theory, the method for direct approximate solution of differential equations, and the 
iteration method. 9 ' 10 The results obtained by the above methods were checked numerically. 

(B) Determination of the electron and phonon temperatures and electric fields in the sample boundary 
areas. Here the above-mentioned methods were used and the following problems were studied: 
determination of temperature fields near the interface of two media and heat transfer through thin layers, 
temperature, and electric fields in planar laminated media. ' L 

For cases A and B, the temperature distribution of electrons and phonons was obtained. 



3.3.2 Low Carrier Concentrations 



If the current carrier concentrations are not very high, the electron temperature approximation cannot be 
used; the respective criteria are given in. 8 The kinetic equation fc 
function /in a one-dimensional approximation is of the form 



used; the respective criteria are given in. 8 The kinetic equation for the calculation of electron distribution 



m ax ap x 

where the operator S e combines the interaction of electrons with phonons and point defects. 
Expanding the electron distribution function into the Legendre polynomials : 

f(p) = X/z(e)P,(cos 0) ( 3 - 17 ) 

1=0 

where ■& is the angle between the direction of the electron momentum vector and the axis x. 
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Expanding all the terms in Equation 3.16 into the Legendre polynomials and using recurrent relation- 
ships, the infinite chain of equations which was solved by the iteration method 1 and the variation 
method can be obtained. The following infinitesimal parameters were used for the average lengths: 



3 Ll 



gT=TJT« 1 (3.18) 



Both methods provided a possibility to determine the electron distribution function. 

3.4 Distribution Function 

Both methods for solving kinetic equations enable structurally equivalent equations for the symmetric 
part of the current carrier distribution function under large temperature gradient. In the first order by 
parameters 3.18 these can be represented in the form 1 ' l 



/^ 



"^ + %(s)T t + %(S)T ^ 



/o -/o 1 + *i(e) -r + ^(b)E — + %(s)T— + %(s)T- T (3.19) 



where fg is the equilibrium Maxwell function, and ^(e) are known functions (presented in 
Refs. [1,2,13,14]). The resulting distribution function has two principal differences from the equilibrium 
Maxwell function 1,2,13 ' 14 : 

First, it depends on the squares of the local temperature gradient and the gradient of electrochemical 
potential. Therefore, the distribution function has nonlinear properties. 

Second, it also depends on the coordinate derivatives of the temperature gradients and on the electric 
field. Therefore, the resulting distribution function is of a nonlocal nature; it depends not only on the 
generalized forces at a given sample point, but also on their values at neighboring points. Obviously, the 
nonlocality and the nonlinearity of the distribution function will create a series of new effects originating 
from a sharp spatial inhomogeneity in a temperature field. 



3.5 Temperature Fields and Transport Coefficients 

All the methods discussed result in equivalent expressions for the electron and phonon temperatures 
under given thermal boundary conditions. The basic results can be illustrated as a consequence of a 
mismatch between the electron and phonon temperatures © — T e — T p . It turned out that temperature 
mismatch in the bulk of a sample can be represented as a sum of three terms © = @ F + © r + ©j, where 
the term @ F appears in the absence of electric current and is related only to the processes of thermal 
conduction (the authors call it "the Fourier mismatch"), @ r occurs at the emergence of heat flow and 
electric current ("the Thomson mismatch"), and ©j exists even without a heat flow, but in the presence of 
electric current (we shall call it "the Joule mismatch"). The analytical form of the temperature mismatch 
terms and their evaluation can be found in Refs. [1,2,10]. 

The temperature mismatch in the volume and at the interface of a sample leads to new phenomena 
which are thermal analogues of the Joule and Peltier effects (Table 3.2). 



TABLE 3.2 Analc 


igues of the Joule and Peltier Effects 








Electric Current (Density;) 




Heat Flow (Density q) 




Sample Volume 


Contract of Two Samples 


Samp 


e Volume 


Contract of Two Samples 


Joule effect 


Peltier effect 


Analogue 


of Joule effect 


Analogue 


of Peltier effect 
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Using the distribution function (3.19), analytical expressions for current and heat flow densities at 
different scattering mechanisms can be obtained. The relations for the flow densities completely 
determine all transport phenomena in the electron subsystem, taking into account the nonlinear and 
nonlocal effects, caused by large temperature gradients. Designating any kinetic coefficient by K, and 
including under K the resistivity p, the Peltier coefficient is 77, the thermoelectric coefficient a, the 
electron heat conductivity K e , and also the average electron energy e. A change in kinetic coefficient K due 
to a large temperature gradient can be represented in the form 1 ' 2 



8Ki = ^Y^ = ^( M «G 2 + u^Gpj + u^p 2 ) 2 ) (3.20) 



where K is the kinetic coefficient under small temperature gradient, index t identifies the method of 
solving (f = p, v, tw for the perturbation theory, variation, and two-temperature methods, accordingly), 
and G ~ (T 2 — T{)l2a is the average temperature gradient. The values of function U, can be found in 
Refs. [1,2]; this function has the order of the small parameter (Equation 3.18) U, ~ g T . The values of 
numerical coefficients u^ are presented in Table 3.3. 

The table gives the values u£ for the cases when nonlocal effects are absent (*c p = const.) and for the 
most important case of temperature dependence (k p ~ T^ 1 ). It can be seen that all the three methods 
give equal signs of corresponding parameters u t '^, and the variation and two-temperature approaches 
lead to a full coincidence of numerical values of these parameters. 

In contrast to corresponding coefficients of standard theory, the kinetic coefficients in the case being 
considered depend on the combinations of temperature gradient and current density squares, i.e., they 
are nonlinear. Second, they are related to high derivatives of temperature with respect to coordinates, and 
hence they are nonlocal. The nonlinearity and nonlocality of the kinetic coefficients result in a series of 
deviations from the standard relationships. Thus, the Lorentz factor L = K e p/T is no longer a constant, 
but becomes a function of current density and coordinate derivatives of temperature. Besides, in the case 
under consideration the thermoelectric Thomson relation 77 = aT is not satisfied. The relationship is 
simultaneously one of the Onsager relations which are valid when there is a linear relation between the 
flows and the generalized forces. Therefore, under large temperature gradients the Onsager relations are 
not valid either. 



TABLE 3.3 Kinetic Parameters Obtained by Various Methods 

Physical Method of u { , 1} 
Parameter Solving (f) 



k~ — const K n 



e (average electrons energy) p —0.83 0.17 —0.33 —0.33 0.67 

0.5 -0.75 -0.75 0.5 

0.5 -0.75 -0.75 0.5 

p (resistivity) p -0.21 0.26 1.04 1.04 0.61 

0.25 -0.38 -0.38 0.25 

-0.37 -0.37 -0.37 0.75 

II (Peltier coefficient) p -1.43/a 2.07/a -0.271a -0.27/ a 21a 

6.25/a -51a -51a II a 

0.5 1.5/a + 0.5 -l/a-0.75 -21a - 0.75 0.75/a + 1 

a (thermo-EMF) p 0.97/a -0.01/a 



0.75/a 
K e (heat conductivity p — 0.76 — 0.26 

of electrons) v — 7 — 7 

0.25 



p 


-0.83 


V 


-0.5 


tw 


-0.5 


p 


-0.21 


V 


-0.25 


tw 


-1.87 


P 


- 1.43/ff 


V 


4.25/a 


tw 


- 1.5/a - 


p 


0.97/a 


V 





tw 


- 0.75/a 


P 


-0.76 


V 


- 7 


tw 


0.25 
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3.6 New Effects 



The abnormal thermoelectric effect of EMF generation (or the Benedicks effect) was investigated. ' 

The Peltier coefficient under large temperature gradient includes the terms which are linear and 
quadratic with respect to temperature gradient. Thus, one can speak about the inverse Peltier effect 1 ' when 
current density changes its direction into the opposite, with invariable temperature field. This is a "pure" 
effect caused by a large temperature gradient. Because large temperature gradients reduce the symmetry of 
crystals, tens of new anisotropic thermoelectric, thermomagnetic, and thermal effects come into 

U • 17-21 

being. 

As long as large temperature gradients change all material properties, the material's requirements for 
thermoelements operating under large temperature gradients also change. ' It turned out that under 
large temperature gradients the role of thermoelectric figure-of- merit is performed by a new parameter 
which is a function of temperature gradient. 

The heat up of carriers under the influence of large temperature gradients brings about a set of new 
transport effects: three mechanisms of electron and phonon temperatures mismatch; the heat analogues 
of the Joule and Peltier effects; the abnormal and inverse thermoelectric effects; and a number of 
anisotropic effects, etc. The nonlinear and nonlocal effects which have been discussed are obviously 
manifested in thermoelectric structures with submicron and nanosize inhomogeneities, such as 
quantum-dimensional objects and thermoelements with point contacts. It is obvious that large 
temperature gradients are realized in these cases, and a rigorous theory of thermoelectric phenomena 
should be constructed with regard for the nonlinearity and nonlocality. It is probable that these large 
temperature gradients are one of the reasons for the reported considerable increase in thermoelectric 
figure-of-merit. 

3.7 The Outlook for Practical Applications 

The experimental creation of large temperature gradients with a view to check new effects is rather compli- 
cated. A relatively large temperature gradient is formed in a "pin-plane" contacts area (Figure 3.2). 
This model was used for the experimental investigation of two pure effects caused by large temperature 
gradients. The abnormal thermoelectric effect of EMF generation (or the Benedicks effect) ' and the 
inverse thermoelectric effect are generally known. The presence of the inverse thermoelectric effect 




FIGURE 3.2 Scheme of a large temperature gradient formation. 1, 2, thermostats with temperatures T t and T 2 , 
respectively; 3, flat homogeneous sample cut along its boundaries 4. 



i 2006 by Taylor & Francis Group, LLC 



Thermoelectric Phenomena under Large Temperature Gradients 



3-9 



h 



1 vr 2 vr. 



2. 



7 



- i i - 



FIGURE 3.3 Scheme of thermoelement with a large temperature gradient and a system of heaters and coolers: 
1, periodic system of heaters; 2, periodic system of coolers; 3, thermoelement working medium; R, external load 
resistance; 2i, heaters temperature; T 2 , coolers temperature; AT 1; AT 2 , temperature gradients. 



corresponds to transition from the case when temperature gradient and current density in the 
semiconductor are parallel to the case when they are antiparallel. Such transition causes a change in 
electric resistivity: 



pij, dT/dx) ^ p(-j, dT/dx) 



(3.21) 



Both effects have been observed in n-Ge, the theory being in fairly good agreement with 
experiment. 15 ' 16 

Thermoelectric devices utilizing the Benedicks effect are of interest. " Figure 3.3 shows a linear variant 
of this thermoelement. A periodic system of heaters and coolers, arranged symmetrically, creates in a 
semiconductor plate periodic areas with large and small temperature gradients. These thermal conditions 
provide a series connection of the Benedicks EMR The system heaters and coolers can move relative to 
each other. This causes a change in the Benedicks EMF direction and, accordingly, creates another 
functional dependence of the thermoelement EMF on the motion. 

Figure 3.4 is a schematic of a spiral thermoelement utilizing a large temperature gradient. The 
thermoelement comprises heaters and coolers located on a spiral in pairs. The EMF across the spiral ends 
is proportional to the number of its coils, the number of pairs of heaters and coolers, the value of the 
Benedicks effect, which in turn depends on the properties of thermoelectric material, and the magnitude 




FIGURE 3.4 Scheme of a spiral thermoelement with a large temperature gradient: 1, spiral working medium; 
2, cooler; 3, heater; T lt heaters temperature; T 2 , coolers temperature; \T X , AT 2 , temperature gradients; R, external 
electric load. 



i 2006 by Taylor & Francis Group, LLC 



3-10 



Thermoelectrics Handbook: Macro to Nano 




7j ^ 1 



(b) 




FIGURE 3.5 Models of thermoelements with a large temperature gradient: (1) large temperature gradient region; 
(a) configuration of module leg with a neck; (b) leg with pressure points; (c) thermocouple with a connecting plate 
(2) forming a large temperature gradient in contact regions; (d) rotating discs with different temperatures. 

of temperature gradient. Because of the absence of connections, such thermoelements are particularly 
stable in time and can be used in high-precision measuring techniques. 

Large temperature gradients can also contribute considerably to a thermoelectric figure-of-merit. 
A value of ZT = 1.7 at room temperature was obtained in thermoelements with pin-plane type point 
contacts. 4 

A series of patents 5 describe the structures of thermocouple devices used to achieve a large temperature 
gradient (Figure 3.5). In some of them temperature gradients up to 4.2 X 10 8 K/cm were achieved. The 
author of the patents claims that the use of such devices allows efficiency values close to the Carnot cycle 
efficiency to be achieved. 
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